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the dynamics of a classical particle moving in the presence of the Thomas-Fermi potential. The 
dynamics is regular for atoms and some dots, but in general in the case of dots, the motion contains 
a chaotic component. We compute the correlation effects. They appear at the order NlnN for 
atoms, in agreement with available data. For dots, they appear at the order N. 
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I. INTRODUCTION 



Ever since the invention by Thomas [1] and Fermi [2] of a simplified theory of an atom 
with a large number N of electrons, many efforts have been made to systematically improve 
on it [3]. Corrections were made, then refined, by Dirac [4], Scott [5], Schwinger [6], and 
Englert and Schwinger [7], which add terms of order TV 6 / 3 and N 5 ^ 3 to the N 7 ^ 3 Thomas- 
Fermi (TF) energy. TF theory and its corrections, collectively referred to as the statistical 
atom, thus seemed to result in an expansion of the atomic ground state energy in iV -1 / 3 . 
After Lieb and Simon [8] proved that TF theory becomes exact when N — > oo, this expansion 
was put on a rigorous basis by Fefferman and Seco [9], who proved that the energy of the 
statistical atom and the energy of Hartree-Fock (HF) theory are equivalent and exact up to 
order iV 5//3 . A crucial step further was made by Schwinger and Englert [10], who showed 
that there exist oscillating corrections to the ground state energy of order N A I 3 and period 
iV 1 / 3 . They interpreted such corrections as indicating shell effects. 

It remains to determine, however, at which order in N both the statistical atom and HF 
theory break down, and to compute the dominant correlation effects, which are ignored in 
both approaches. It is this task that we undertake in this paper. 

But we are also interested by the same type of problem in the case of quantum dots. We 
consider in this paper quantum dots to be 2-dimensional artificial atoms, whose N electrons 
are subject to a confining potential and interact by way of the standard 3-dimensional 
Coulomb interaction. The determination of the ground state energy as a function of N is 
of great interest especially since, for a class of confining potentials, the energy has become 
experimentally accessible [11]. The analogy between atoms and quantum dots became quite 
clear when Lieb, Solovej, and Yngvason [12] proved that in the case of dots also, TF theory 
becomes exact when N — > oo. We are therefore interested in this case also by the corrections 
to TF theory. An important difference with the atomic case is that the confining potential 
isn't necessarily rotationally symmetric. This leads to shell effects sensitive to the nature of 
the dynamics of a classical particle moving in the presence of the TF potential. Depending 
on the chosen confining potential, this dynamics can be fully regular, fully chaotic, or, most 
frequently, mixed. Therefore, quantum dots are ideal laboratories for the study of quantum 
chaos. 

To determine the ground state energy of large atoms and dots, our main idea comes 
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from the fact that after a simple Z-dependent rescaling of the coordinates — where Z is the 
number of protons in the case of atoms, and N = Z in the case of dots — the Coulomb 
interaction between the electrons becomes a weak and long range interaction, indicating the 
validity of a mean field theory when N, Z — > oo [13]. But at the same time, the role of 
the Planck constant is played by a parameter e = Z^ 1 ^, d being the dimension of space. 
Therefore, large values of Z correspond also to a semiclassical limit. In our case, the mean 
field theory is simply Hartree theory, and its semiclassical limit is TF theory, in agreement 
with the theorem of Lieb and Simon [8]. 

In order to go beyond TF theory, our strategy is the following. Considering first e as 
an independent parameter, we derive an asymptotic expansion in Z~ x of the energy. The 
dominant term of this expansion is Hartree theory, and the corrections to it correspond to 
exchange and correlation effects. But remembering that in our case, e is small, we then 
take the semiclassical limit of each term in this expansion. In this way, we have computed 
the correlation contributions to dominant order. In the case of atoms, they give a term 
aZ In Z + bZ , and in the case of dots a term cZ. The other contributions to the energy 
correspond to a Hartree-exchange (HX) theory, which coincides with HF theory up to a 
certain order in Z (at least Z 5 ^ 3 in the case of atoms). 

In the case of neutral atoms (Z — N), we compare our results for the correlation energy 
with experimental and numerical values. It was suggested [14, 15] on the basis of these 
values that the dominant term of the correlation energy is of the order iV 4 / 3 , contrary to our 
results. But we can see that our logarithmic correction allows to fit well the data. 

In the case of dots, we completely determine the smooth part of the HX energy to the 
order N. The oscillating part is less important than in the case of atoms, and its analysis 
is deferred to another article. In the case of atoms, after summarizing known results, which 
come from the HX energy, we indicate what remains to computed for this part of the energy. 
This is a delicate problem in semiclassical physics that we have not solved. Our results show 
that the Schwinger-Englert oscillations can be interpreted as resulting from a trace formula 
for an integrable system. The integrability in this case is due to the rotational symmetry of 
the TF potential. 
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II. THE MODEL FOR ATOMS AND QUANTUM DOTS 



We use the dimensionless Hamiltonian 

N N 1 N 1 

ff=-^A J + ^%.) + -^-— (1) 

3=1 3=1 3^k \ X 3 Xk \ 

to describe an atom or quantum dot containing N electrons. In the case of atoms, we've 
used half the Bohr radius for the unit of length and two hartrees for the unit of energy; a 
similar choice can be made for dots. In the case of atoms, Z denotes the number of protons 
and V(x) = — \x\~ 1 . We have neglected relativistic effects. In the case of dots, Z will be 
identified with N, and V(x) is a confining potential whose form is not specified. At this 
stage, V(x) is independent of N, but if we describe a specific experiment, we may have to 
consider a smooth dependence of V{x) on N. 

After the rescaling of the coordinates x = Z 2 / d ~ 1 x, where d is the space dimension, with 
d = 3 for atoms and d = 2 for dots, the Hamiltonian H becomes H = Z 2 ~ 2 ^ d H, where the 
new Hamiltonian H is given by 

N N N 

^E^+Evm+^E^- (2) 

j= i j= i 1 J fc| 

Accordingly, the ground state energies E of H and E of H are related by E = Z 2 ~ 2 / d E. In 
H, e = Z^ 1 ^ plays the role of H, and the Coulomb interaction between the electrons looks 
like a mean field type interaction when N — > oo, considering that linijv->oo Z/N — 1. Taking 
first e to be an independent parameter, we will give an exact formula for the ground state 
energy E as an expansion in the small parameter Z~ l . This formula decomposes E into two 
parts, the first corresponding to the HX energy, the second to the correlation energy. We 
then evaluate the two parts for small e, determining first the leading order of the correlation 
energy, then considering the expansion in e of the HX energy up to this leading order. 

III. GRAND-CANONICAL FORMULATION 

To determine the ground state energy E, we start with the grand-canonical partition 
function Q((j,,(3), with \i the chemical potential, and (5 the inverse temperature. From 
Q(fj,,(3), we can get the ground state pressure P(n) as 

P(fu)=\im ilnQ^/J). (3) 



If we define by 

E» = fiN-P(fi), (4) 
then the ground state energy at a given number N of electrons is given by 

E = va£E> t . (5) 

In the case of atoms, the presence of a continuous spectrum of H would give an infinite 
Q(fl, //), but the problem is easily solved by adding to V(x) a confining potential, suppressed 
after the limit /? — ^ is taken. For clarity of notation, we don't explicitly write this additional 
confining potential. 

Using coherent state representation, we can write the grand-canonical partition function 
as [16] 

Q = J Z#*,V] exp-^y dxJ^dT^{x,T) [d T -fi + h']^(x,r) 

ip a (x, r), ip*{x, t) being Grassmann variables and h! the one-body Hamiltonian hi = — e 2 A + 
V. By applying a Hubbard- Stratonovich transformation on the Coulomb interaction in this 
integral, we can integrate over the Grassmann variables, so that the partition function 
becomes 

Q = (Det 2 [d T -n + h' + iZ' 1 ' 2 ^] ) _^ (7) 

(f>(x,r) being a Gaussian field of zero mean and covariance C~ 1 (x,r x \y,T y ), C being the 
operator of kernel C(x, r x \y, r y ) = 5{t x — r y ) \ x — y| -1 . Here, and in what follows, we denote 
the determinant and trace Det and Tr if they operate on both space and time, and det and 
tr if they operate on space only. 
Let us now make the shift 

<f>{x,T) = -iZ 1/2 w(x) +9(x,t). (8) 
The partition function can then be written as 

Q = Q (exp [iZ 1/2 (C^w, 6)] Det 2 [l + iZ 1/2 K6] ) , (9) 



d;C 

where 

Q = Det 2 [d T -n+ h] exp \\z(w, C^w)] , (10) 





h being the new one-body Hamiltonian h = h' + w, and K the operator K = [d T — fi + h] 1 . 
On the other hand, 

Det 2 [1 + iZ-^KB] = exp -2 £ Tr (^)" (11) 

n=l 

We can find a w(x) such that 

• z i/2 (c-i lu> e) + 2iZ~ 1 / 2 Tr = 0, (12) 
so that the linear term in 6 disappears. Indeed, 

/ (13) 



w (X) 



where 



Po(x) = —n(x\x), (14) 



with 

n(x\y) = (x, i + ^ {ph) y), (15) 

does the job, since —K(x, r\x, r) = n(x\x). 
If we introduce the operator T of kernel 

T(x,T x \y,T y ) = 2K(x,T x \y,T y )K(y,T y \x,T x ), (16) 

we see that 2Tr(K9) 2 = (0,r#), and we can write the partition function in the form 



where 



Q = Q Qi(expA(6)} , (17) 



^("gC'^ (18) 



and 



-l) 



nZ n/2 
n=3 



9 being a new Gaussian field of zero mean and covariance C~ x — T/Z in (17) and C _1 in 
(18). But Q\ is simply 

Q 1 = [ Det (1 - CY/Z)] - 1/2 = exp - £ — Tr(CT)™. (20) 

n=l 
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We therefore decompose the pressure P(/x) into four terms as 

P = P H + P x + P C;1 + P C;2 , (21) 

where 

P H = lim -mQ , (22) 

is the pressure in a mean field approximation, w(x) being the mean field. The remaining 
terms 

P x = lim ^Tr(CT), (23) 

1 °° 1 

P c .! = lim — V Tr(CT) n , (24) 

n=2 

and 

P c . 2 = lim -m/exp,4(0)> . (25) 

correspond to fluctuation effects around the mean field. This decomposition of the pressure 
will correspond to a natural one for the ground state energy, and the indices H, X, and C 
foreshadow the nature of the corresponding contributions to the energy. 



IV. GROUND STATE ENERGY: HARTREE-EXCHANGE AND CORRELATION 
DECOMPOSITION 

In correspondence with the decomposition (21) of the pressure, we decompose the ground 
state energy as 

E" = E* + E* + E^+E^. (26) 

.Eh is given by 

E£ = pN - P H = pN - 2tr [{p - h)#(p - h)} - f / ^L po ( x ) Po (y), (27) 

^ J \ x y\ 

where we have used equation (13) which defines w(x), and where po(x) is given by (14), 
but with both w(x) and po(x) now obtained from the density matrix n(x\y) in the ground 
state n(x\y) = (x\'&(p — h)\y). We can write E£ in the form 

E£ = pN-2 [" l deD(e) - § / ^Lp ( x ) Po (y), (28) 
J 2 J \x - y\ 



7 



D(e) being the integrated density of states D(e) = tr$(e — h) of the Hamiltonian h, which 
we now write as h = — e 2 A + W(x), where 

W(x) = V(x) + w(x) = V(x) + [ -^— Po (y) (29) 

i r y\ 



1 f dxdy 



is the self-consistent potential. 
Ex is given by 

E " = " Px = ~Z,I ]x^~y~\ nl ^ yl (30) 
and describes exchange effects. E^,. x and E%. 2 describe correlation effects. For E^ ;1 = — P C;1 , 
it is useful to introduce the representation of the kernel of K 

K(x,r x \y,r y ) = l-J2 K ^ x \y) eMTx ~ Ty) > ( 31 ) 

where 

K u = (iu-fi + h)- 1 , (32) 

being the Matsubara frequencies uj = n{2n + for n G Z. Then if r x 7^ r y , we have 

the representation of the kernel of Y 

T(x,r x \y,r y ) = - £ r n (a ! |y) e <n < T "- T »\ (33) 

with 

r n (x|i/) = i^iir w (x|i/)iir w+ n(i/|x). (34) 

In this way, we see that 

^Tr(CTr = i^tr(CT^, (35) 

where in the right hand side of this equation, the kernels of the operators C and Tq, and 
the trace are defined on M. d . The limit (3 — > oo is simply taken by replacing -| ^ n in (35) 
and ^ XL m (34) by J ^ and J respectively. Consequently, 

^ = -iEi,/i'OT. (36, 

n=2 

Finally, 

= -^c ;2 = Km - j ln(exp A(0) >, ;C _ 1 _ r/z . (37) 

Up to this point the expressions given for the different terms of the energy compose 
an exact asymptotic expansion in Z^ 1 for E^ . From now on, we'll consider only those 



terms which will contribute, in the subsequent semiclassical limit e — > 0, to the dominant 
correlation energy. For this purpose, whereas we'll need to keep all terms in E£ ;1 because of 
the subtlety of its semiclassical limit, we'll only need to keep E%. 2 up to order Z~ 3 (with e 
fixed). If a n = -2(-i) n /n and A n (9) = tr(K9) n , then 



n=3 



and we write 



-E%. 2 = lim — 



0,4 



2^3\ 3/e ; c-i T zsV"°/0;C- 



1^ 



where 



5 = i™ Z [(^> e . C -i_ r/z " <^4> 0;C _ 1 ]. 



(38) 



0(Z- 4 ), (39) 



(40) 



V. FIXING THE CHEMICAL POTENTIAL 

There exists a yet unknown parameter in the expressions of the energies, namely the 
chemical potential ji. We therefore need to determine the dependence of \i on N. We can 
write E^ in the form 



1 °° 1 



n=2 



where 



and 



P Q (fi) = 2 J\eD(e) + ^ J ^^p ( x )p (y), 



Pi 



<"> = / 



dxdy 9 



F - y| 

the terms -P n (/-0 describing correlation effects when n > 2. 
The chemical potential is fixed by the equation 



OO 



n=0 



Writing the chemical potential in a Z expansion as 



Z™' 

n=0 



(41) 



(42) 



(43) 



(44) 



(45) 
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and the energy as 

oo 

E = V 

Z 



n=0 



and keeping only the terms necessary to determine the energy up to the leading 
correlation energy, we get 

E = p N-P (p ), 

E 1 = fjuN - HiPo(no) - PiM, 

E 2 = p 2 N - frP'oM - l^P"M ~ PiP'M - P2M. 
If we fix p by the constraint 



we see that 



and therefore 



PjM 



E = p N - P (p ), 
E 1 = -PM, 

E p(u) 

Let us note that equation (50) can be written in the more natural form 

N = 2D(p ). 

Indeed, 



d 

dp 

and 



^ deD(e) = D(p) — tr [$(p - h)d„W] , 



/dxdy f dxdy 
rPo(x)p (y) = 2 -p Q (x)d fl p (y), 
\x-y\ J \x-y\ 

but from equation (29), 

j \ x y\ 

so that P^p) = 2D(p). 

We see now that the energy E H , defined by 



E K = E» 

n=no 
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corresponds to the Hartree energy, because equation (55) fixing /i corresponds to the nor- 
malization condition 

r n 

(60) 



/ 



dxp (x) = ^, 



of the electronic density po(x) associated to h. 
We can now write the ground state energy as 



E = E HX + E c , 



(61) 



where E HX is the Hartree-exchange (HX) energy 



E HX = E H + E x + A, 



(62) 



with the exchange energy 



and 



E X = E£ 



A = I w 



Ml 2 

x 



and where E c is the correlation energy E c = E c -i + E Cl2 , with 



and 



E c -1 = £c:l 



-^c;2 - ^c;2 



M=M0 



(63) 



(64) 



(65) 



(66) 



As we will now only need the dominant term /i of the chemical potential, we will denote in 
what follows /i = /i . 



VI. THOMAS-FERMI THEORY 



The TF energy corresponds to the dominant contribution to the ground state energy. It 
is obtained by keeping the dominant term in Z -1 , which is Hartree theory, then taking the 
semiclassical limit e — > 0. As shown in (28), the Hartree energy is completely determined 
by the knowledge of the integrated density of states D(e) and the density po(x). The 
semiclassical limit of D(e) is given by 

D(e) = ^ J dx$(e- W(x)), (67) 
11 



and, introducing the local chemical potential 



fi(x) = fi-W(x), (68) 

and the notation 

fx + (x)=^x)$(fi(x)), (69) 
the semiclassical limit of po(x) is given by 

p (x) = 2a d J dx^Tix), (70) 

where a d = Sd/(dxd), with Xd = (2vr) d and Sd the surface of the <i-dimensional unit sphere 
S d = 2n d l 2 /T(d/2). The normalization condition in the Hartree approximation (60) becomes 

2a d J dx/j, d J 2 (x) = —, (71) 

which results in the self-consistent equation for p(x) 

fi(x) = p- V(x) - 2a d J ^y—^f{y). (72) 

Combining these results, we get for the semiclassical limit of E H the TF result, expressed 
for the original Hamiltonian H, 



= a d Z 3 -^ 2 



fi + J dxn d J 2 (x)V(x) + J dxfi^^ix) 



(73) 



In the case of atoms, fi(x) = p(\x\) because of the rotational symmetry of V(x), and we will 
use the fact that fJ>+(x) = 0(\x\~ r ) when \x\ — > and [J,+(x) = (9(|a;| _4 ) when \x\ — > oo. 
And we can recall that fi = in the case of neutral atoms (Z = N). In the case of dots of 
confinement such that V(x) — > oo when \x\ — > oo, it has been proven [12] that /i(x) has 
compact support and is bounded on this support. 



VII. SEMICLASSICAL ESTIMATE OF THE CORRELATION ENERGY 

We now consider the correlation energy in the semiclassical limit. In particular, this limit 
corresponds to taking, in E c -i and E C;2 , for the chemical potential \x and the self-consistent 
potential W(x), their TF values, as defined in the previous sections. 
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A. Energy E C]1 



From (32) it follows that 



K u {x\y) = de- 



Pe{x\y) 



too — ji + e 



(74) 



where p e {x\y) = (x\5(e — h)\y). Therefore, since 



duj 



T n (x\y) = 2 / —K u (x\y)K u+n (x\y), 



we have 



T n (x\y) = / de 1 de 2 



I 

p ei ( x \y)pe2(y\x) 



ifl + e\ — e 2 
From (36), we see that E C]1 is given by 



i3(e 2 - p) - 7?(ei - p) 



where 



oo 

2 71 
n=2 



Let us first consider the dominant term in Z 1 , namely Y" 2 , and define E'*^ 
can write as 



F 2 



d\id\<2 



tr(G Al G 



A2 / 



+ Ai)(ifi + A 2 )' 



(75) 
(76) 

(77) 

(78) 
-F 2 /4. We 

(79) 



where 



G\(x\y) = 2 / rfeirfe 2 / rfz C(x\z)p ei (z\y)p e2 (z\y) 

x 5(A - (ei - e 2 )) f#(e 2 - /i) - (d - //)] . (80) 



Then 



z 2 



or, more conveniently, 



Ai — A 2 
2 /,0 ° 

Z 2 " 



/•oo 

/ dt tr(L 2 ), 
Jo 



with 



L t (a;|i/) = 2 ^ d ei de 2 ^ dzp ei (z| 2 /)p e2 (z| 2 /)C(a;|z)e- t ( ei - e2 ^(e 1 - p)#(p - e 2 ). 



(81) 



(82) 



(83) 
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Inserting the semiclassical expression for p e (z\y) 

p e (x\y) = ±- ! d p5{e-W{^)-p 2 )< 

€ Xd J 



,ip(x-y)/e 



(84) 



we get 



L t (x + er/2\x-er/2) = j dqdk I dz C(p - k/e)e ixq - izq a t (k] /i(z)) , (85) 

^ Xd J J 



where C (p) is the Fourier transform of the Coulomb potential 



C(p) = ^\p 

Xd 



and 



at(k;fi) = J dq 1 dq 2 e- t ^-^$(qj- l 2)$( l 2-q 2 2 )5(k + q 1 + q 2 ). 
In this way we obtain the useful representation 



o2 

d-1 °d j£ 



where 



and 



ix\ e ' 

X e = J dpdkkA(k,p)\k + ep/2| 1_d |fe - ep/2] 1 - 



A(k,p) = 



dt 



dx a t (k; fi(x))e 



ipx 



Xdh Jo 

where the k~ x factor is chosen so that A is finite when k — * 0. 



(86) 



(87) 



(88) 



(89) 



(90) 



Let us try to apply the semiclassical limit simply by setting e = in X e . We have 

4 r°° r 

dpA(k,p) = -J dt J dx a 2 (k, p(x)) (91) 



and since a t satisfies the scaling relation a^fc/i 1 / 2 , /i) = p d / 2 a^ t (k, 1), we find, using equation 
(71), 

(92) 



A — — 77 1 d: 



S, Z 



where 



dt J dkk 2 - 2d a 2 (k; 1). 



(93) 



In the case of dots, I 2 is finite and, after a somewhat lengthy computation [17], we can 
evaluate it to 

J 2 = 27r 3 (l-hi2), (94) 
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so that with N = Z, 



= -(l-hi2) = 0.1534. 



(95) 



It is interesting to note that this result is universal, that is, independent of the confining 
potential. 

In the case of atoms, corresponds to the first second order contribution, of r s per- 
turbation theory, to the correlation energy for the homogeneous electron gas. It is loga- 
rithmically divergent. For the homogeneous electron gas, a finite result of order lnr s for 
the correlation energy is obtained by summing the most divergent terms of all higher order 
contributions [18]. The atom, however, isn't fully homogeneous, it is effectively confined. 
This confinement, related to an inhomogeneous chemical potential, translates into an effec- 
tive cut-off in the integral in (89), which will result in a finite X € of order lne -1 . A lengthy 
computation (see appendix) gives the result 



X e =X log lne- i +X lin + (9(e), 



(96) 



where X log = 2(1 - ln2)(27r) 5 iV/Z, and 
X lin = ^(2tt) 5 - ^ + | In 2 - 2G(1 - In 2) - 4 In 2 2 



i(47r) 4 [(l-ln2)A-i-fi . (97) 



where G is Catalan's constant, 



poo 

A= dr r 2 fi 3 / 2 (r) In fi 1 / 2 (r) , 
Jo 



and 



with 



B 



poo poo 

= dt dp ln(p)g 2 (p,t), 
Jo Jo 



g(p,t) = / drrfi' + (r)e~ tfl + (r) cos(pr) — - 

Jo L pr 



(98) 



(99) 



(100) 



Note that the logarithmic term is universal, but the linear correction contains non-universal 
terms. 

It remains, however, to consider the remainder 



2^n 



(101) 
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in the semiclassical limit. It is convenient to rewrite R as 

1 ^ (-l) n f dVL 



and 



n=3 

where 

X e n (Q) = tr [ - Z^C 1 / 2 ^ 1 / 2 ] (103) 

because it easily follows from (76) that the operator —Z~ 1 C 1 ^ 2 T 6 qC 1 ^ 2 is self- adjoint and 
positive definite. We can prove that in the case of dots, lim e ^ e f dQX^(Q) = 0, for n > 3, 
so that E c -i = E*. t . But in the case of atoms, J d£lX®(£l) is finite when n > 3. 

We now consider exclusively the case of atoms. We can put the semiclassical representa- 
tion of the kernel of -Z^C 1 ' 2 ! ' en C 1/2 in the form 

[-Z- l C l / 2 Y en C l ' 2 ]{x\y) = -i— / dk^T^k^e^-^, (104) 

where 

T^k 1 \k 2 ) = ^Jdye^ k ^%(k 1 + k 2]f , + (y)), (105) 

bh(P^) = -J dq n2 + ll] q y H(tP/4 + q) 2 - ^ - (ep/4 - qf) ■ (106) 
Therefore 

*n(«) = ^ tr (?S) B - (107) 

This representation is useful because fr^(p; /x), given by 

b°n(P^) = J dq6(<?-ti n fo£* qr (108) 

is finite. More explicitely 6^(p; /i) = /i 1//2 6(p/i 1 / 2 /fi), where 6(x) = 27r[l — a;" 1 arctanrr]. 
It remains to verify that X®(Q) is integrable. For this purpose, we will use the inequality 

Km < [x° 2 m n/2 - But 

X 2 °(fi) = 2 f d Zl dz 2 [ y ^e^-^n^Cg;^^-)). (109) 

J J [P + Q) [P- Q) jJi 

Consider first the case VI 3> 1. Writing X$(Q) as 

X2 ^ ) = lj ip+qyZ-q)* I ^ e ^?(P^ + ( a; +^)) & ?^^( £C -^))' ( 11Q ) 
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we see that 



lim QX^Q) = 2x1 I dx^ 2 (x) 



l Jo P 



dp h2{ \ 

2 (?) 



(111) 



which is finite. Hence X^fl) = 0(Q 1 ). In the limit Q <C 1, we start from the expression 

X 2 °(n)=45f r%G Q {p). (112) 
Jo P 



Here 



and 



Gn(p) 



dx g(x) I dr pfi 1 / 2 (r)b 



1/2/ \ 

(r, 



r sin(xpr) 



,,11, x+1 
g{x) = ; rln 



One can show that 



and 



x x 2 + 1 | a; — 1| 
Jo P 3 

fpM<c 2 , 
^ dp 



(113) 

(114) 

(115) 
(116) 



(117) 



Jn P 3 

if Q < Qq. Consequently, if Q > Q , < di/ft"/ 2 , and if Q < Q , \X°(Q)\ < 

d 2 {lnQ- 1 ) n / 2 , so that / §X°{tt) is finite if n > 3. 

We can now summarize all these results. In the case of atoms, we have 



N 



E C]1 = 0.03109 ^ In Z 1/3 + Z~ 1/3 [0.03700 + x a + x b + x] , 



(118) 



where x a = 0.01979 A, A being defined by (98), x b = 0.01027 5, B being defined by (99), 
and 



x 



4e 



2 " n 

n=3 



2vr trl ^ j U(2tt)5 



;ii9) 



Tq being the compact operator of kernel 



T°( kl \k 2 ) = ±-J cte^f(*) e ^-^) 6 (J*i±M^/ 2 (z)), (120) 



where &(x) = 2n[l — x 1 arctanx]. 
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B. Energy E c - 2 

Starting from (39) and anticipating the result of the semiclassical limit, we decompose 
E c - 2 as 

9 

-E c , 2 = 2a 4 A d + a 4 B d + -a\C d + R d . (121) 

In order to give a simplified expression for these terms, we introduce the notations (xj, Tj) = 
j, C(x j ,r j \x k ,r k ) = C jk , and K(x j ,r j \x k ,r k ) = K jk , and define L 12 ... n = Yl] =1 with 
K n (n+i) = K nl . Then 

i r 4 

A d = lim — T[dj L 1234 C 12 C 34 , (122) 

/3->oo JZ Z J - LJ - 

.7 = 1 



J 

4 



and 



where 



B d = lim -J- [f[dj L 1234 C 13 C 24 , (123) 

C d = lim -L / dld2/!/ 2 , (124) 
h= [ d2d3L 123 C 23 . (125) 



If 



then 



and 



R d = \a\D d + a A E d + a 6 F d , (126) 



1 /" 6 

D d = lim — - / TT dj L 123 L 456 C 14 \c 25 C 36 + C 26 C 35 ] , (127) 

/3^oo pZ° J ^ L J 



E d = lim -L- / <ftd2(CTC) 12 AK 12 (KKK) 21 + 2 / dSdA K 13 C 34 K 41 K 32 K 24 , (128) 
where = KjkCjk, and 

^ ^ j=i fc=i 

The semiclassical computation of R d , which corresponds to terms of order Z~ 3 requires a 
lot of work, that we do not reproduce here. The final result is that R d = 0(e 2d ~ 3 lne" 1 ), so 
that we can ignore R d , since the other terms will be of order e d ~ 2 . 
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We begin by first giving the expressions for Ad, Bd, and Cd when the limit (3 — > oo is 
taken. They are 

Ad = ~^J ^Z^^i^Pe^Pe,)^- ^(^-e 2 ), (130) 
where now K(x\y) = — n(x\y)C(x\y), 

B d =4~2 L i 6l,e3;e2 ' e4) ^-^(ea-^-e^-e^, (131) 

Z 7 -M- ei + e 3 - e 2 - e 4 



where 



L(ei,e 3 ;e 2 ,e 4 ) = ^ Y\_ dx j C ( x ^-\ x ^) C ( x 2\ x d 



and 



where 



X pe 1 ( X l\ X 2)pe i ( X 3\X4)pe 2 (X2\x 3 )p e4 (x 4 \x 1 ), (132) 



C ' = ^/ra^W- (133) 



/(x) = -2 J ^^(p„kp a )(x\x)#( ei (134) 
Let us then take the semiclassical limit of these expressions. We have 

-K( Xl \x 2 ) = -rl— [ dqe^- x ^g d {q,p{^)), (135) 

e Xd J 



where 



9 d (q,p) = c d J ^(^"(q-p) 2 ), (136) 

with c d = Sd/Xd, so that 
A d 2 



Xd 



dq 1 dq 2 dp l dp 2 dxdsdrdtg d {qi + f ; p(x + e s ^))g d (q 1 - f; p(x + e^)) 



x e *[-(q2,«)+(p2,r)+(t,qL-p2)] ^((P^ + f ) 2 - + ffM^j - (Pi ~ f )^ , (137) 

2(9, P2) + - M 33 + es ) 
It is now easy to take the limit e — > 0, so that 

lim-^ = -— / ^p^(p;//(*))5(p 2 -//(*))■ (138) 
The scaling relation p~ 1 ^ 2 g d (p 1 ^ 2 p; p) = gd(p; 1) and the normalisation condition (71) give 
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since g d (l; 1) = A/S d . 

To compute B d , we write L(ei, e 3 ; e 2 , e 4 ) semiclassically as 



L(ei, e 3 ; e 2 , e 4 ) 



x gi(T'i,pi+q2)+i(T , 2,P2-2qi) 
, 2 



dp 1 dp 2 dq 1 dq 2 I dxdridr 2 \pi + ^l 1 ^Ipi-^l 1 d 



x % - (qi + f) - W(a: + ^))5(e 3 - (<h - |) - W(a: - e*)) 
x <5(e 2 - (qr 2 + f ) 2 - W(x + ef ))<J(e 4 - (q 2 - f ) 2 - W(x + ef )) . (140) 

In this form, the limit e — > is easily taken, and we get 
,2 d c 2 



B d = e 



d-2- 



* / ^%iPi+p 2 r-iPi-p 2 r- d / dx/dh.^M*)), (i4i) 

Xd J Pi~P 2 J 



where 



f( Pl ,p 2 ;fi) = J dq#((q-p 1 ) 2 -t,)#((q + p 1 ) 2 -t,) 

X^-(q + p 2 ) 2 )^-(q-p 2 ) 2 ). (142) 



But /i- d / 2 /(/i 1/2 Pi,/i 1/2 p 2 ; 1) = /(Pi,P 2 ; 1) so that finally 



B d = e' 



2 



where the constant b d is given by 

dpidp 2 



b d = 



v\-vl 



Pi + P2MP1 - p 2 | 1 - d /(Pi,p 2 ;i). 



(143) 



(144) 



Remarkably, this integral appears in the second-order exchange contribution to the corre- 
lation energy of the homogeneous electron gas [18]. First computed numerically, its value 
was then obtained in closed form for d = 3 by Onsager et. al. [19], for d — 2 by Isihara and 
Ioriatti [20], and for any d by Glasser [21]. Effectively, 



B* = 



eN 



5 ln2 -i« 3 » 



and 



B 2 = ^L 



3 7T 2 J ' 



where G is Catalan's constant and 



" ^ (n + l) 3 ^ 2m + 1* 



(145) 



(146) 



(147) 



ra=0 



m=0 
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To compute Cd, we write f(x) semiclassically as 

- f( x ) = f dk 1 dk 2 dqdzdrg(q; fi(x + ez )) e i(™-fci)-i(*,fc2) 

e X'd J 

^((fci + | ) 2 ~ + e(f + + 6(f ~ |)j - (fci - | ) 2 ) 

2(fc 1 ,fc 2 )+^(a J + c (f-5))-|i(x + c(f + =)) • 1 } 

The limit e — > in this expression gives 

-f( x ) = ^r~ I dqg(q;fi(x))5(q 2 -^x)), (149) 

e Xd J 

so that 

Q = e- 2 4 / ^^rf- 1)/2 (a ;i )/,r i)/2 (- 2 ). (150) 
X d J \xi-x 2 \ 

It is a priori surprising that such a term of order Z~ 3 gives semiclassically a result of the same 
order as those of order Z~ 3 . Indeed, all the other terms of order Z~ 3 gave semiclassically a 
result of the order Z~ 3 e d+3 , up to lne^ 1 corrections. We assume that all the term of higher 
order in Z^ 1 give semiclassically a result of the order Z~ n e d+n , up to corrections in lne -1 . 
But we have not proven it, and this question remains to be settled. 

To summarize, in the case of atoms, and using the rotational symmetry of TF local 
chemical potential, we find 

N 1 2f27r) 3 f°° f r 
-E c , 2 = 0.06390 ^ - zl/3 Y / drr J dss 2 fi + (r)fi + (s), (151) 

and in the case of dots, we find 

-E c , 2 = 0.1455 - -L f ^^Vl\ Xl) ^\ X2 ). (152) 

Z7T J \X\ — X\ 

VIII. CORRELATION ENERGY OF NEUTRAL ATOMS: COMPARISON WITH 
NUMERICAL AND EXPERIMENTAL VALUES 

Let us first recall that we have decomposed the ground state energy as E = E HX + E c . 
While E HX looks like the HF energy E HF , the two energies differ in their inclusion of exchange 
effects, and shouldn't be confused. Taking a determinant of Hartree wave functions for the 
trial wave function, we will get E HX for the energy. Therefore, we have the inequality 
-^hf < £ H x- As shown by a semiclassical HF theory, the difference E HX — E HF is zero up 
to the order iV 5 / 3 [14]. We have therefore shown that HF is correct up to order iV 5 / 3 , a 
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0.07 r 




Figure 1. Per electron HF-relative correlation energy for neutral atoms with up to 55 electrons. 
Data points correspond to experimental values (exp) [22] and extended HF values (ext HF) [15] for 
E£ F /N. For sufficiently large atoms, the per electron energy E* F /N = -0.062 In TV - 0.018 (solid 
line) agrees with both experimental and extended HF values, presenting an essentially oscillating 
deviation of less than 8% when N > 10. 

result that has already been rigorously proven by Fefferman and Seco [9]. However, we do 
not know whether E HX and E HF are equivalent up to order N. Therefore, we do not know 
if E c corresponds to what is commonly referred to as the correlation energy, that is the 
HF-relative energy E^ F = E — E HF , up to order N. 

Furthermore, while we have numerically computed the terms x a and xj, of E C; i as well as 
the second term of the contribution E c - 2 in (151) for neutral atoms — resulting in contribu- 
tions to -E c , given in hartrees, of 0.06533 N, -0.00329 N, and -1.1044 N, respectively- 
we haven't computed the constant x appearing in the contribution E c .^. Therefore, for neu- 
tral atoms, we have the correlation energy, given in hartrees, E c = —0.062 N In N + cN, with 
c to be determined. If we now assume that E HF differs from E HX at most by a contribution 
5N, we can compare the energy E^ F = —0.062 N\nN + c'N, with d = c + 5, to experi- 
mental and numerical values for E^ F . Experimental values exist for atoms containing up to 
18 electrons [22], and were obtained by removing from measured ground state energies the 
relativistic contribution; numerical values exist for atoms containing up to 55 electrons [15], 
and were obtained in an extended HF approach. We see in figure 1, that the NlnN term in 
E c is essential in reproducing the behavior of reference values, and that with d = —0.018, 
Eq F agrees with experimental and numerical values, differing essentially by an oscillating 
contribution when N is sufficiently large. 
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IX. HARTREE-EXCHANGE ENERGY OF QUANTUM DOTS 



Our remaining task is to compute the HX energy semiclassically, in the case of dots. It 
has been recognized [23] that one should distinguish a smooth part and an oscillating part 
in a semiclassical expansion of the density matrix or the integrated density of states of a 
quantum system. The oscillating part, contrary to the smooth part, depends crucially on 
the nature of the classical dynamics associated to the potential. We will consider here only 
the smooth part, the oscillating one being discussed in another article. 

The density matrix has the semiclassical expansion 

n(x + er/2\x - er/2) = e~ 2 n a (r; p(x)) + n b (r; p(x)) + 0(e), (153) 

where 

n a (r;fi) = J dpe^(p-p 2 ) (154) 

For n b (r,p), we will only need here the fact that 

J dxn b (0;p(x)) = J dx AW(x)S(p(x)). (155) 

Consequently, we have, since D(p) = J dxn(x\x), 

J'deDie) = dxplix) - J dx AW{x)d{p{x)). (156) 

This suggests the decomposition of p and W, up to the order e 2 , 

p = p a + e 2 /if>, (157) 
W(x) = W a (x) + e 2 W b (x), (158) 



and correspondingly 



with 



and 



po{x) = po-,a(x) + e 2 p -b{x), (159) 
po-a( x ) = 2n a (0;p a (x)), (160) 

p b (x), (161) 



p - b (x) = 2n b (0; p a (x)) + 2-^-n a (0; p) 



where p a (x) = p a — W a (x) and p b (x) = p b — W b (x). p a will be fixed by the constraint 



1 = J dxp ., a (x), (162) 
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and /jLj, by the constraint 

= J dxp - b (x). 
The self-consistent equation decomposes into 

w a (x) = v(x) + [ y^—pUy), 
J \ x — y\ 

and 



W b (x) 



dy 



iPo-M- 



(163) 



(164) 



(165) 



\x - y\ 

We can now decompose the Hartree part E H of the ground state energy using the fact that 
Z = N 

E H = E H . a + e 2 E H . b , (166) 



where 



N f N f 

E H , a = p a N- — J dxp 2 a (x)$(p a (x)) -- J dxp . a (x)(W a (x)-V(x)), (167) 



and 



E S;b = N 



/i b - dx p . a (x)p b (x) - / dx p Q . a (x)W b (x 



N 



j dxAW a (x)^ a (x)). 



N f 

+ — J dxAW a (x)&(fi a (x)) 



24tt 

The exchange term becomes 



(168) 



L f^L f 



dx 



n a (r; p a {x)) + e 2 n b (r; p a {x)) 



or 



where 



/:, -J^ I dxp 3 Jl(x) + 0(e), 



J = 



e (2tt) 3 
dp±dp 2 



Pi +P2\ 



16tt 



Finally, the correction 



A= i_[d,E x \ 



2 P»(jm) 



(169) 



;ito) 



(171) 



(172) 
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gives 



A = 



J dx^J 2 (x)a{x) 



3 r ' 
71 / dx a(x)-&(ijL a (x)) 

where a(x) is the solution of the linear integral equation 



(173) 



(174) 



We can now summarize the results for the smooth part of the ground state energy corre- 
sponding to HX. Expressed for the original problem, we have 



E BX = N 2 E TF + N 3 ^ 2 ^ J rf^f(x) 



N f N 
+ — I dxAW(x)$(fi(x)) + 



24tt 



dx fj 1 / 2 (x)a(x) 



n J dx a(x)'d[fi(x)) 



-, (175) 



where E TF is the TF energy 



E — - 



7^ J dxV(x)fi+(x)], 



(176) 



and ji and W(x) being solution of the TF equations 



W(x) = V(x) + — 



dy 



and 



2-7T J \x — y 



1 ■= — / dxfi+(x). 



(177) 



178) 



X. HARTREE-EXCHANGE ENERGY OF ATOMS 



In the case of atoms, it is needed to evaluate semiclassically the HX energy to the same 
order as the correlation energy. One can take advantage of the spherical symmetry of the 
potential W(x), thus reducing the problem to a one-dimensional one, which can be studied 
semiclassically by WKB type techniques. However, a specific difficulty arises in the atomic 
case. The potential W(r) diverges like the Coulomb potential r~ l near the origin. In physical 
terms, the semiclassical limit has to be reconsidered for strongly bound electrons. This is 
at the origin of the Scott and Schwinger corrections, which have successfully handled this 
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problem. But one needs more, namely to compute the integrated density of states up to 
the order e. Moreover, the exchange energy, whose dominant term has been computed by 
Dirac [4] requires a knowledge of the density matrix n(x\y) up to the order e 2 . If we use 
a standard expression for this correction of the density matrix, we get a logarithmically 
divergent correction for the exchange energy. The origin of this divergence is the slow decay 
of the Coulomb potential appearing in the expression of the exchange energy. 

As in the case of dots, the energy will be decomposed into an oscillating and a smooth 
part, E HX = E HX . S + E HX . osc . We expect for E HX . S a neutral atom an asymptotic expansion 
given by 

7 

E HX . s = J2 c i Nj + c ° NlnN > ( 179 ) 

i=3 

with the constants, expressed in hartrees, 

c 7 = -0.7687, c 6 = -0.5, c 5 = -0.2699, (180) 

known [3]. We have not undertaken the task of computing C4, C3, and c . 

The oscillating part, of order iV 4//3 , has been computed by Schwinger and Englert [10]. 
In the case of atoms we see therefore, that contrary to the case of dots, it is more important 
than the correlation energy. This oscillating part can be understood as the first appearance 
of shell effects, in an atom described as a liquid by the smooth part of the energy. This 
is the interpretation of Schwinger and Englert. But it also has a dynamical interpretation. 
Indeed it is standard now to decompose the density of states d(e) = d e D(e) into two parts 
semiclassically [23] as 

d(e) = d s {e) + d osc {e). (181) 

The smooth part d s (e) is given by an asymptotic expansion in e _1 , whose coefficients are 
some integrals depending on W(x). The corresponding part of E HX was given in (179). The 
oscillating part 

rfosc(e) = 7) cos (-S(e, 7) + <r 7 |) (182) 

7 

is given by a sum over the periodic orbits 7 of a classical particle moving in the potential 
W(x), where S(e,^) is the classical action along the orbit, 07 is the orbit's Maslov index, 
and A e (e, 7) depends on the orbit's period and stability. d osc (e) therefore depends crucially 
on the nature of the dynamics associated to W(x). In the case of atoms, the TF potential 
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W(x) is rotationally symmetric, so that the dynamics is that of an integrable system. In 
this case the periodic orbits are stable and degenerate. A general formula for rf OS c(e), in the 
case of a rotationally symmetric potential, is given in [24]. One can see that it corresponds 
to the formula obtained by Schwinger and Englert. In fact, all these authors base their 
analysis on the use of the Poisson formula 



XI. CONCLUSION 

We have provided a method to systematically compute the ground state energy of non- 
relativistic atoms and quantum dots, by means of an asymptotic expansion in N^ 1 , whose 
terms have have to be evaluated semiclassically. The dominant terms are give by a semi- 
classical HX theory, which coincides with HF theory only up to a certain order in N' 1 . 
Correlation effects go beyond HF. They are of order N In N for atoms, and order N for dots. 
It remains in the case of atoms to compute numerically a constant appearing in the term 
of order N in the correlation energy. In the case of atoms, it remains also to fully compute 
the HX energy up to the order N. This represents a challenge in semiclassical physics. This 
computation would allow a better comparison with the data, because when we used our 
results for the correlation energy, this quantity was defined as the difference between the 
true energy and the HF energy, and we have seen that it is not the full HF energy which 
matters, but rather the HX energy. 

Among the possible extensions of this work in the case of atoms are the following: (a) For 
non-relativistic atoms, and with L the total angular momentum and S the spin, L 2 , L z , S 2 , 
and S z are conserved. It would therefore be interesting to compute the ground state energy 
with these quantities being fixed; (b) Compute the ground state energy taking into account 
the dominant relativistic corrections, which become more important when N is large. 

Experimental results for quantum dots in the presence of a magnetic field have been 
obtained [25]. A TF type theory has already been established in this case [12]. It is therefore 
an interesting, though challenging, problem to extend our results for the corrections to TF 
theory to this situation. 




(183) 
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XIII. APPENDIX 

In the case of atoms, the computation of E* yl is rather delicate, so that we give here a 
summary of the main steps in this computation. The crucial quantity to compute is what 
we called X e , and given by 



with 



//»oo 
dkdp (k + ef )' 2 (k - e?y 2 J dt \T t (k\p)\' 

T t (k\p) = J dxe ipx a t (k;fi(x)), 



where we recall that 



a t (k; fi) = J dq e-^ q+k)2 - q2] ${(k + qf - fj)^(fi - q 2 ) . 



r t (k\p) depends only on k and p, so that 

X e = M J dpY £ (p), 

where 



(P)= / 
Jo 



k 2 

dk— In 



k 2 + l 



k + 1 



k-1 



7(^|P), 



where 



7(fc|p) J dt . 



We then decompose Y t into four parts. Let 



dk 



k 2 



k 2 + 1 



In 



k + 1 



k-l 



- d(k - !)■ 



7(ef|p) 



Y, A p) = 2ff 1 (e k i\p) 

7 (e^|p)- 7 (0|p) 
W)=2f f 7 (0|p). 

J evil ft 



Y 3 ,(p) = 2 1 — 



1 dk 

ep/2 



(184) 
(185) 

(186) 

(187) 

(188) 
(189) 

(190) 
(191) 
(192) 
(193) 
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This decomposition is justified by the fact that it can be shown that y(0|p) is finite and 
7(fc|p) = 0(k~ 3 ). For these proofs one needs to remember that fi(r) ~ r~ l when r — > and 
/i(r) ~ r -4 when r — > oo. The decomposition allows us to get the asymptotic behavior of 



X 



where 



and £> = J^ =1 Bj, with 



and 



X e =Alne _1 + S, 



A = 16n J dp7(0|p), 



C 



£1 



8ttC / rfp7(0|p) 
"°° dk 



16tt 



dp J y7(^|p) 
167T | rfp^ 1 ^[ 7 (A;|p)- 7 (0|p)] 
4tt y dpln(2/p) 7 (0|p), 



A: 2 



A; 2 + 1 



In 



A;+ 1 



k- 1 



(194) 
(195) 

(196) 
(197) 
(198) 
(199) 

(200) 



other work is needed, that we do not reproduce here but can be found in [17], to express 
these constants in terms of /i(r), then obtain expression (97) by using the normalization 
condition (71) where applicable. 



[1] L. H. Thomas, Proc. Camb. Phil. Soc. 23, 542 (1927). 
[2] E. Fermi, Rend. Accad. Naz. Lincei 6, 602 (1927). 

[3] J. D. Morgan, in Atomic, Molecular, & Optical Physics Handbook, edited by G. W. F. Drake 

(AIP Press, 1996) Chap. 10, pp. 233-242. 
[4] P. A. M. Dirac, Proc. Camb. Philos. Soc. 26, 376 (1930). 
[5] J. M. C. Scott, Philo. Mag. 43, 859 (1952). 
[6] J. Schwinger, Phys. Rev. A 22, 1827 (1980); 24, 2353 (1981). 

[7] B.-G. Englert and J. Schwinger, Phys. Rev. A 26, 2322 (1982); 29, 2331 (1984); 29, 2339 
(1984); 29, 2353 (1984). 

29 



[8] E. H. Lieb and B. Simon, Phys. Rev. Lett. 31, 681 (1973); Adv. in Math. 23, 22 (1977). 
[9] C. L. Fefferman and L. A. Seco, Adv. in Math. 107, 1 (1994); Helv. Phys. Act. 71, 1 (1997). 
[10] B.-G. Englert and J. Schwinger, Phys. Rev. A 32, 47 (1985). 

[11] L. P. Kouwenhoven, D. G. Austing, and S. Tarucha, Rep. Prog. Phys. 64, 701 (2001); S. M. 

Reimann and M. Manninen, Rev. Mod. Phys. 74, 1283 (2002). 
[12] E. H. Lieb, J. P. Solovej, and J. Yngvason, Phys. Rev. B 51, 10646 (1995). 
[13] J. L. Lebowitz and O. Penrose, J. Math. Phys. 7, 98 (1966); E. H. Lieb, 7, 1016 (1966). 
[14] B.-G. Englert, Semiclassical Theory of Atoms, Lecture Notes in Physics (Springer- Verlag, 

1988). 

[15] E. Clementi and G. Corongiu, Int. J. Quantum Chem. 62, 571 (1997). 

[16] J. W. Negele and H. Orland, Quantum Many-Particle Systems, Frontiers in Physics (Addison- 

Wesley Publishing Company, 1988). 
[17] R. Rueedi, Ground State Properties of Large Atoms and Quantum Dots, Ph.D. thesis, Ecole 

Polytechnique Federale de Lausanne (2009). 
[18] M. Gell-Mann and K. A. Brueckner, Phys. Rev. 106, 364 (1957). 
[19] L. Onsager, L. Mittag, and M. J. Stephen, Annalen der Physik 473, 71 (1966). 
[20] A. Isihara and L. Ioriatti, Phys. Rev. B 22, 214 (1980). 
[21] M. L. Glasser, J. Comput. Appl. Math. 10, 293 (1984). 

[22] S. J. Chakravorty and E. R. Davidson, J. Phys. Chem. 100, 6167 (1996). The values for the 
HF-relative correlation energy are located in the supplemental pages, available online 

[23] M. Brack and R. K. Bhaduri, Semiclassical Physics, Frontiers in Physics (Addison- Wesley 
Publishing Company, Inc., 1997). 

[24] M. V. Berry and K. E. Mount, Rep. Prog. Phys. 35, 315 (1972). 

[25] S. Tarucha, D. G. Austing, T. Honda, R. J. van der Hage, and L. P. Kouwenhoven, Phys. 
Rev. Lett. 77, 3613 (1996). 



30 



